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ABSTRACT 

We have measured spatial and temporal variability in the y band sky brightness over the course 
of four nights above Cerro Tololo near Cerro Pachon, Chile, the planned site for the Large Synoptic 
Survey Telescope (LSST). Our wide-angle camera lens provided a 41 ° field of view and a 145" pixel 
scale. We minimized potential system throughput differences by deploying a deep depletion CCD 
and a filter that matches the proposed LSST band (970 to 1030 nm). Images of the sky exhibited 
coherent wave structure, attributable to atmospheric gravity waves at 90 km altitude, creating 3% 
to 4% rms spatial sky flux variability on scales of about 2 degrees and larger. Over the course of a 
full night the j/3 band additionally showed highly coherent temporal variability of up to a factor of 
2 in flux. We estimate the mean absolute sky level to be approximately 2/3 — 17.8 mag(Vega), or 
1/3 = 18.3 mag(AB). While our observations were made through a y^ filter, the relative sky brightness 
variability should hold for all proposed y bands, whereas the absolute levels should more strongly 
depend on spectral response. The spatial variability presents a challenge to wide-field cameras that 
require illumination correction strategies that make use of stacked sky flats. The temporal variability 
may warrant an adaptive y band imaging strategy for LSST, to take advantage of times when the sky 
is darkest. 

Subject headings: atmospheric effects, instrumentation: photometers, techniques: photometric, spec- 
troscopic 



1. INTRODUCTION 

The Large Synoptic Survey Telescope (LSST) 
is a proposed wide-field camera that is cur- 
rently in the design and development p hase 
(ILSST Science Collaborations: Paul A. Abell et all 
120091 hereafter LSST09). The project has designated 
the telescope's site to be Cerro Pachon, Chile, on the 
same ridge that hosts the Gemini South and SOAR 
telescopes. The LSST will survey the sky in the Sloan 
Digital Sky Survey's ugriz passbands, plus an additional 
near infrared band, y, at wavelengths ~ 990 nm, only re- 
cently made possible by deep depletion CCD technology 
that enhance s quantum efficiency (QE) at the reddest 
wavelength s (10 'C onnor eTall [20061 ) . The Pan-STARRS 
(jKaiser et al. 2002) survey also uses the y band. 

The designation of "y" may present confusion. In- 
frared astronomers use a band F, typically with a cen- 
tral wavelength near 1.03/^m HHillenbrand et all 120021: 
[Hewet t et al. 2006). Their use of infrared detectors, with 
a flatter response curve near a micron, yields a passband 
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with significantly different shape than the CCD-based y. 
Our filter also bears no relation to the optical y filter of 
the Stromgren passbands, which is centered on 550 nm. 
This degeneracy in terminology is most unfortunate, but 
seems likely to persist. 

The sky brightness in the y band has not been well 
characterized for astronomical CCD applications. Past 
work describing sky brightness at observato ries typically 
cite values for the UBVRI bands, but not ?/ (jPatat II2003I : 
IKrisciunas et al.ll2007tlSanchez et al.ll2007t) . The reddest 
band in the SDSS survey is z, whose red edge is deter- 
mined by the rapidly falling QE of the SDSS detectors, 
which have very little sensitivity in the y band region. 
The UKI RT survey presented sky measurements for the 
Y band (jWarren et al.l l2007f ) , but as explained before, 
the total transmission curve is significantly different than 
that using a deep depletion CCD. 

Our aim has been to provide the first dedicated CCD 
measurements, at the proposed LSST site, of the y band 
sky background and its variability over the course of a 
night. While the exact character of the new y band to 
be used for LSST is under active consideration, we used 
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Fig. 1. — The LSST y band and its variants (top), and the 
sky emission spectrum (bottom). The top panel shows the 
expected transmission spectrum for proposed y band variants 
ULSST S cience Collaborations: Paul A. Abell et al.l 1200 91 includ- 
ing predicted optics and detector response, as well as the relative 
transmission that we measure for our apparatus (black solid line). 
Our effective filter curve closely resembles the planned yi band. All 
spectra exclude atmospheric absorption and scattering, which we 
plot separately as a dashed line. The bottom panel shows a typical 
sky spectrum in photons m~^ nm"'^ a rcsec~'^, dominated by 
Meinel band OH emission linos (taken from lPuxlev et ani200lf ). 

a filter candidate called (LSST09). We expect our 
qualitative findings to hold for any y band variant un- 
der active consideration because the sky background in 
all cases arises overwhelmingly from narrow OH emis- 
sion lines. Figure [T] shows potential y band transmission 
curves with detector quantum efficiency (QE) and atmo- 
spheric response for both LSST and our apparatus, along 
with the anticipated sky emission spectrum. 

In !j2]we briefly describe the known nature of sky emis- 
sion within the y passband, covering past aeronomic 
( §2.ip and astronomical ( §2.2p measurements. In [J3]we 
describe the dedicated instrument we assembled for these 
measurements. Section 2] discusses the observations, and 
the data reduction is outlined in fj5] Section |6] presents 
results on spatial (i j6.1[) and temporal ( ii|6.2l) variability, 
and an absolute calibration (i i6.3|) . and we conclude with 

m 

2. THE NATURE OF Y BAND SKY BACKGROUND 

The background sky in the y passband is domi- 
nated by Meinel band emission due to hydroxyl (jMeinell 
in a physical proce s s that has long been 
observed (eg 'Wi ens et al.l 11 997: Nakamvu'a ct al. 1999") 
and modeled (iRousselot et al. .2000.; Cosby et al, 2006; 
iCosbv fc Slangeii I2007D bv atmospher ic scientists (for a 
useful overview, see lMarsh et al.ir2006[ ). Hydroxyl is well 
strati fied in the mesopause at 90 km in a 10 km thick 
layer CB aker fc Stahlll988l : iMelo et al.ll2000( ). Above the 
ozone layer in the stratosphere, daytime UV light from 
the sun dissociates O2 to make 0-1-0, which react with 
other O2 to produce O3. UV light also dissociates O3, 
completing the oxygen-ozone cycle. O3 additionally re- 
acts with H to produce O2 and excited OH, which then 
emits photons to produce the airglow within the y band. 



As distinct from Rayleigh scattering, Meinel emission 
is known to vary significantl y over a number of spa- 
tial and temporal scales (eg iMarsh "erall[2006l) . The 
airglow varies with the tides, and is linked to the sea- 
sons and the solar cycle. Traveling and standing atmo- 
spheric gravity waves, which are coherent density pertur- 
bations analogous to ocean waves, are regularly observed 
by aeronomers ( eglTaylor fc Hil]||T99lH G" arcia et alSlWA 
iNakamura et al.lll999l: iLi et al.ll2005l) . and give rise to a 
variable background over a single, wide-field y band CCD 
image. 

2.1. Aeronomic Measurements 

Aeronomers, like astronomers, have made extensive use 
of CCDs since their invention, but also have not ex- 
plored the y passband region with them in depth be- 
cause the development of deep depletion is recent. CCDs 
have been used to measure sky emission lines in isola- 
tion with Fabry-Perot interferometers, for the purpose of 
monitoring atmospheric gravity waves and variability in 
species temperature and emission ra te (eg, the Spectra l 
Airglow Temperature Imager, SATI. lWiens et al.lll9"97[ ). 
They have also been used to image gra vity waves di- 
rectly with narrow and broad fi lters (eg iTavlor fc Hilll 
[l99lt iGarcia et al.|[l997l: ILi et al.ll2005l) . These applica- 
tions have been most often used in concert with radar, 
lidar, balloons, rocketry, and space-based experiments to 
achieve fully three-dimensional dynamical information of 
O2 and OH atmospheric layers. None of these observa- 
tions ventured far into the red of today's deep depletion 
CCDs, nor did the analyses achieve calibrated photome- 
try. Existing aeronomic data are therefore of limited use 
to LSST, as LSST requires prior knowledge not only of 
the variability but also of the expected absolute back- 
ground level in the y channel. Our instrument and data 
reduction techniques aim to provide all needed informa- 
tion at once. 

2.2. Previous Astronomical Measurements 

Previous astronomical work is also insufficient, as the 
y band sky has not been directly characterized. Extrap- 
olating to the y band using other astronomical filters 
is problematic because of the highly structured nature 
of the Meinel band emission (Figure [l] bottom panel). 
Small changes in the filter's red edge may suddenly alter 
the bac kground level. For example, the recent UKIRT 
survey (jWarren et al.ll2007l ) does include a Y band, but 
using an infrared detector and filter whose total re- 
sponse function do not match that expected for LSST. 
iWarren et al.l (|2007D in fact showed that their Z band 
sky brightness from the UKIRT observations on Mauna 
Kea was 0.4 mag arcsec^^ brighter than that obtained 
by SDSS in their z band, at Apache Point, New Mexico. 
Besides differences in effective filtering, the discrepancy 
may be explained by differences in atmospheric condi- 
tions in time and in space given the great geographical 
separation. We therefore must assume that the UKIRT 
median sky brightness oiY = 17 A mag arcsec"^ (Vega) 
cannot be applied to LSST. 

Our effective filter function closely matches LSST's t/3 
band and thus minimizes potential errors due to differ- 
ences in relative spectral response; we additionally ad- 
dress the geographical problem by observing at Cerro 
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Tololo, a mountain ridge 10 km northwest of LS ST's 
planned site, Cerro Pachon. iKrisciunas et al.l (|2007[ ) pre- 
sented a historical overview of sky brightness observa- 
tions from Cerro Tololo, for the period 1992-2007. They 
claimed that there was no discernable difference in sky 
brightness between that site and Cerro Pachon. This 
leads us to conclude that our y band sky brightness mea- 
surements from Tololo should be applicable to the LSST 
site on Pachon. 

We conclude that previous work within aeronomy and 
astronomy is not sufficient for characterizing the y band 
sky background for LSST, and that the precise nature of 
the chosen y filter affects the observed background due to 
the multitude of strong OH lines at ~ 1 /Ltm and redder. 
The dependence of the apparent NIR sky brightness on 
the filter and detector combination, and the potential for 
variation between observatory sites, motivates the mea- 
surements described here. 

3. APPARATUS 

We used a Pixis system (detector, readout electron- 
ics, dewar and shutter) from Princeton Instruments, 
equipped with a type 1024BR CCD. This detector, 
made by EEV with their designation 47-10, was a back- 
illuminated deep depletion sensor that exhibits enhanced 
sensitivity at wavelengths near 900 nm due to the thicker 
region of silicon from which photocharge is harvested. 
The sensor was an array of 1040 x 1027 pixels, each of 
which was 13.0 /im on a side. The Pixis camera was 
equipped with a thermoelectric cooler that maintained 
the detector at a temperature of —20 C. Dark current 
was negligible compared to the flux seen from the sky. 
We measured the camera's gain to be 3.8 ADU^^. 

Princeton Instruments designed the 1024BR CCD to 
suppress fringing that would otherwise arise from night 
sky line emission. We observed no significant fixed fring- 
ing patterns in any of the images we obtained, and did 
not have to apply any fringe-frame corrections to the 
data. We therefore judge any systematics in the photom- 
etry due to fringing to be subdominant for the purposes 
of the results presented here. 

The Pixis camera was controlled from a laptop com- 
puter that ran the Linux operating system. We used 
custom software to control the exposure times and imag- 
ing sequence of the instrument. All images were stored 
as 16 bit FITS files. 

We attached a 17-55 mm variable focal length f/2.8 
Canon EPS lens to the Pixis dewar, using a C-mount 
adapter made by Birger Engineering. The system was 
operated at f/2.8. The y filter we used was an interfer- 
ence filter manufactured to the LSST 7/2 specification by 
Barr Associates, in a 2 in x 2 in format. We attached the 
filter to the front of the lens using a modified professional 
photography filter holder that mated to the 75 mm front 
threads of the Canon lens. The lens's focal length was 
adjusted so as to provide a plate scale of 145" pixel" ^, 
and a square field of view of 41°. Although we make 
use of a 2/2 filter, the resultant system throughput more 
closely matches the LSST 2/3 response curve due to the 
difference in our detector QE and antireflection coatings 
on the optics from the predicted LSST system. 

The camera, lens and filter combination was mounted 
on a stationary tripod. This arrangement did not pro- 
vide pointing information in the image headers, but as 



described below we were later able to determine astrom- 
etry using the positions of cataloged sources. 

Our optical configuration was not identical to that 
of LSST. The LSST fihers will be placed in an //1.23 
converging beam. Light from a point source therefore 
passes through the LSST filter in a (hollow) cone with 
an opening half-angle between 14.2° and 23.6°. Obscu- 
ration from the secondary blocks rays at angles less than 
14.2°. The transmission curve through the interference 
filter depends on the angle of incidence. Rays that tra- 
verse the filter at angles other than normal incidence 
produce a blue shifted transmission function, shifted to 
the blue by roughly 



SX = X 
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where n is the effective index of refraction of the filter 
dielectric. The effective LSST passband is the appropri- 
ately weighted integral of these angle-dependent trans- 
mission functions. 

For our observations we placed the filter in front of 
the lens, so rays from any point source pass through the 
filter as parallel rays. The center of our field is observed 
through the filter at normal incidence. The edges of our 
field are observed through the filter at an incidence angle 
of 20°. While we can't replicate the distribution of rays 
that will pass through the LSST filters, by judiciously 
choosing the region in our field where we measure the 
sky brightness, we can try to match the median LSST ray 
that traverses the filter at an angle of 18.9° (LSST09). 
We therefore expect, even for a perfect match between 
our interference filter and that of LSST, slight differences 
in the effective passbands. See Figure[l]for a comparison 
of passbands. 

4. OBSERVATIONS 

We obtained images over four nights, UT 2007 Sep 13, 
14, 18, and 19. On the 13th we acquired data on the 
second half of the night, and pointed the camera on the 
local meridian and low to the horizon at a zenith angle 
of about 60°. For the remaining nights we aimed the 
camera on the local meridian toward the equatorial, at 
approximately 30° north. 

The moon was 6% illuminated on the first two nights 
and set at 0l''43™ UT. The moon was 44% illuminated 
on the last two nights and rose at 05^30™ UT. For the 
all observations, the moon was not a complicating factor. 
We estimate that the conditions on the first two nights 
were photometric, whereas high clouds began to form on 
the second half of the 18th and all of the 19th. Table □ 
summarizes the observations. 

On each night we elected to take all images at a fixed 
azimuth and elevation in order to simplify our analysis of 
the temporal variation in sky brightness. The nominal 
pointing changed somewhat from night to night, how- 
ever, because we disassembled the apparatus each morn- 
ing. We did not track the camera in R.A. We alternated 
between sequences of five 20 s ("short") exposures and 
ten 300 s ("long") exposures throughout the night. On 
the 14th we acquired a total of 55 short exposures and 
110 long exposures. The short frames had only slight 
trailing of stars, and were useful for astrometry and pho- 
tometric normalization using the normal toolkit of as- 
tronomical data reduction for point sources. The longer 
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TABLE 1 
Summary of the Data. 



Night 


Num 


exposures 


Conditions'^ 




5 min 


20 s 




2007 09 13 


54 


31 


Photometrie 


2007 09 14 


111 


58 


Photometric 


2007 09 18 


76 


40 


First half of night photometric 


2007 09 19 


67 


35 


Clouds 



^ Local date at start of night. 

^ Based on our observations, archived SASCA red-channel an- 
imations (see |http://«M«.ctlo. noao.edu/-sas ca/sasca.htm ), 
GOES East^ 115 micron satelhte imagery (see 
http://www.ctio.noao.edu/sitetests/GQES/), CTIO DIMM 
and Hux monitor data, and Blanco Cosmology Survey observers' 
logs. 

images offered better signal to noise for the sky bright- 
ness, but stars trailed significantly in the EW direction. 

5. DATA REDUCTION 

We used the utilities in IRAF to do a line-by-line bias 
subtraction of each of the images, using the overscan. 
The 300 s exposures were combined into a normalized 
sky flat, by taking a flux-scaled median of the frames. 
This normalized sky flat was then used to flat-field each 
of the images. The normalized sky flat had a fractional 
response suppression of about 20% from the center to the 
corners of the field, presumably due to vignetting and 
perhaps plate scale variation from the lens. This varia- 
tion in fiux sets an upper bound for potential systematic 
errors in sky brightness due to passband dependence on 
the angle of incidence, at 0.2 mag. We take this as an 
extreme upper limit, since we do expect cos" factors to 
roll off the response at the edges of the field. Because 
we are measuring sky brightness, the stacked sky flat is 
indeed appropriate for correcting the pixel-to-pixel vari- 
ations on the scale of a PSF. Our 41° x 41° field of view 
rendered twilight flat-fielding impossible due to sky gra- 
dients over such a large held and to the short window of 
time over which useful sky levels per pixel could be ac- 
quired, and we had no system in place for obtaining the 
equivalent of dome-flats. As noted in fj3l the 1024BR 
CCD was speciflcally designed to suppress fringing, and 
indeed we observed no fringe patterns in any of our data 
and therefore made no fringe-frame corrections. Figures 
[5] and [3] show typical long exposures, after flat-fielding 
and bias subtraction. 

6. RESULTS 

We flrst explored relative measurements of the sky 
background, because we expect relative errors to be sig- 
niflcantly smaller than the overall uncertainty from abso- 
lute zeropoint calibration. We then estimated the overall 
zeropoint in the Vega and AB magnitude systems. 

6.1. Spatial Variability 

The gravity waves we observed in the vast majority 
of our images are easily seen in Figure [5] and [31 These 
structures were not due to flatfielding errors. This can 
be deduced from the two example images themselves, 
which were taken from the same night and had the same 
fiatfield correction applied to them. Roughly speaking, 
the locations of large-scale brightness minima in Figure [2] 
correspond to maxima in Figure [3l Indeed, these waves 




Fig. 2. — A fiat-fielded 300 s exposure in j/3, taken at fixed alti- 
tude and azimuth (righthand frame has stars removed to emphasize 
the airglow structure). The star trails subtend about 1.25° in the 
EW direction; north is up and east is right. This image shows two 
distinct gravity waves traced by NIR narrowband emission, which 
gives rise to short term variability in 5/3 band sky brightness. The 
high frequency ripples have wavelength of about 5 deg and prop- 
agate at between zero and about a degree per minute. The low 
frequency ripple's wavelength is roughly the field of view (40 deg) 
and travels at a degree per minute. Root mean square sky varia- 
tions in a given frame are up to 4% , while the overall sky level over 
a full night varies by a factor of 2. 




Fig. 3. — Same as Figure[2l but taken later in the night, just as 
our calibration star p Ceti entered the field of view (upper right 
corner of lefthand frame). 

appear and disappear, and move coherently in images 
from any given night viewed consecutively. Their direc- 
tion of propagation is perpendicular to their wavefronts. 

To quantitatively estimate the effect of the waves on 
sky brightness at any given time, we computed the mode 
of the sky brightness values for three regions of the detec- 
tor, for each of the flatfielded 300 s images. The reference 
region was the central 500 x 500 pixel area of the detector. 
This was used to estimate the zero-level over the large 
area of sky that we observed. Our "calibrator" region 
was a smaller, circular aperture of radius 50 pixels near 
the detector corner that we ultimately used for zeropoint 
calibration 16. 3p . The "control" region was another 50 
pixel radius aperture near the opposite corner, but at the 
same angle of incidence with respect to the camera bore- 
sight. The calibrator and control regions were chose to 
be smaller than the size of the smallest coherent spatial 
gravity waves in the images, whereas the reference region 
was chosen to be much larger. 

We calculated the mode of pixel values in each aperture 
to attain sky counts c in ADU pixel" ^. Relative magni- 
tudes between each corner region and the central ref- 
erence region were calculated as AS = 2.5 log(ccai/crcf)j 
(likewise for Cctri/crcf ) with units of mag arcsec"^, assum- 
ing a common pixel scale between the two apertures. Fig- 
ure 2] shows these fractional differences over each night. 

Observations on the night of the 13th were taken at 
sec(z) = 2 airmasses at the center of the detector, with 
the calibrator region at 1.4 airmasses and the control re- 
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Fig. 4. — Spatial variation in j/3 band sky brightness versus time, showing 2%-5% rms variability in flux. This plot shows the fractional 
difli'erence in sky brightness between our calibrator region (circular aperture of radius 50 pixels) versus a central reference region of size 
500 X 500 pixels (points), and between a second 50 pixel radius region on the opposite corner of the chip and the reference region (line). 
The large variation on the 13th is due to the extreme difference in airmass between to two corner regions; 5.8 for the northeast corner and 
1.3 for the southwest corner. Early on in the night on the 14th, the Galactic plane fell within the NW aperture, again causing larger than 
typical variation. 



gion at approximately 4 airmasses, due to our large field 
of view. This gave rise to rms spatial sky brightness 
variations of 3% between the calibrator and the refer- 
ence region, 6% between the control and the reference 
region, and 8% between the calibrator and control re- 
gions. On the other three nights, the camera pointed at 
about 1.2 air masses, with the calibrator region at 1.0 air 
masses and the control region at 1.4 air masses. Ignor- 
ing the early part of the 14th, where the Galactic plane 
fell within the control region, these nights exhibited 2%- 
3% rms spatial variability with respect to the reference, 
and 3%~4% rms variability between the calibrator and 
control regions. We conclude that, on angular scales of 
the wavelength of the gravity waves, > 2%, the spatial 
variability in the y band is typically 3%-4% rms. Our rel- 
ative photometric errors are expected to be subdominant 
to this effect. Under the assumptions that gain, exposure 



time, and pixel scale are all constant in any given image, 
these results are independent of those quantities and of 
any absolute flux calibration. 

6.2. Temporal Variability 

The overall temporal sky variability is estimated using 
a similar approach as in the previous section. The fluxes 
within the 500 x 500 pixel reference region are compared 
to the mean over all 4 nights, and plotted in Figure [SJ 
Here we are comparing 2.5 times the log of sky counts in 
ADU pixel" ^ between images rather than within them, 
but because all exposure times, gain values, and pixel 
regions (and therefore pixel scales) were the same, this 
is equivalent to measuring relative sky brightness in mag 
arcsec"'^. Again, all that needs to be done to bring this to 
an absolute scale is to add a zeropoint in mag arcsec"^. 

Immediately evident both in the Figure and in the raw 
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Fig. 5. — Temporal variation in j/3 band sky brightness versus time, showing highly coherent variability of up to a factor of 2 in flux over 
the course of a night. Flux was taken along a fixed line of sight that is 20° from the zenith and directed north. Data from the 13 and 14 
are most probably photometric, the first half of the 18 may be useful, and the 19 is not photometric. The sky brightness levels in this plot 
were obtained from the same 2.13° radius circular region on the detector, and constitutes a differential measurement with very high signal 
to noise ratio. The fiuctuations on both long and short timescales are real. We consider the Vega-based j/3 band zeropoint to be reliable to 
0.2 mag, but this systematic uncertainty is swamped by the temporal variability in sky brightness. To convert to AB mag arcsec"'^, add 

0.5 mag to the 1^2 "^'^ values. 



images viewed in succession is a significant overall sky 
level variability. The peak-to-peak variation on the 13th 
is 1.8 in flux, and on the 14th, 18th, and 19th are 2.0, 2.7, 
and 2.2 in flux, respectively. We have used a huge region 
of about 400 deg^, and we expect our statistical uncer- 
tainties from both the gravity waves and from electron 
counting to be subdominant to the temporal variability 
we measure. Changes in overall sky level are highly cor- 
related in time, and all of the structure seen in Figure [5] 
is due to the atmospheric emission variablity. 

The relative results have made use of only ADU counts, 
assuming constancy for all other quantities that are 
needed for an absolute flux calibration. Our final task 
is to estimate the overall zeropoint. 



6.3. 



measurements from 
^ we determined 



Absolute Sky Brightness 

In order to convert 
ADU s^^ pixel" ^ to mag arcsec 
the gain, pixel scale, and magnitude zeropoint. We 
measured the gain to be 3.8 e" ADU"^ The pixel scale 
and zeropoint measurements were both accomplished 
using a standard star. 

We used photometry of a type AOV standard star to 
establish a conversion between our instrumental magni- 
tudes and a Vega-based ys ba nd magnitude scal e. This 
follo ws the appro a ch us ed in iHillenbrand et al.l (|2002[ ) 
and i Hewett et al.l ()2006D , where AOV stars are used as 
proxies for Vega to establish a zeropoint on the Vega 
magnitude scale. 

We were fortunate to have beta-test access to the 
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TABLE 2 
p Ceti data from the 

LITERATURE. 



TABLE 4 

Systematic error budget. 



Band^ 



^Vcga 



u 


4.85 


0.03 


B 


4.85 


0.03 


V 


4.88 


0.03 


J 


5.08 


0.20 


H 


4.86 


0.02 


K 


4.81 


0.02 



^ U BV da ta arc from 
IMermilliod"! (Il987a) and 
JHK data a re from 
ISkrutskie eFaLl l l20Ml . 



TABLE 3 

p Ceti data. 



Frame 


UT 














(pixel) 


(pixel) 


(ADU) 


(ADU pixel-i) 


110 


06'^19™01= 


925.4 


257.7 


92.6 


14.6 


111 


06'U9™22= 


923.5 


257.8 


117.0 


15.5 


112 


06'U9™44'= 


921.7 


257.9 


129.6 


13.6 


113 


06'^20™05= 


919.1 


257.8 


119.2 


15.3 


114 


06'^ 20™ 26= 


916.8 


258.0 


136.7 


15.3 



^ Centroid of the star. 
^ p Ceti flux. 
Local sky. 



astrometry.net web utility (http://astrometry.net/), 
which has the remarkable ability to determine a WCS 
solution for an arbitrary image. We input one of our im- 
ages, frame 110 (taken at UT 06^19^22'') to this site 
and learned that the image was centered at R.A. = 
0lhl7'"34.8", Decl. = -02°20'49" (J2000). We then 
queried the SIMBAD archive to find AO stars that co- 
incided with this field. The AOV star p Ceti at R.A. = 
02h25"^57.0", Decl. = -12°17'26" (J2000) (also desig- 
nated as HR 708 and SAO 148385) met our criteria. This 
object has cataloged magnitudes that are listed in Table 

H 

We identified p Ceti in five of the short exposure im- 
ages, and used IRAF to perform aperture photometry. 
The resulting centroid positions, object fluxes, and sky 
background values are given in Table [H We deter- 
mined the flux from the star by averaging these mea- 
surements, obtaining 119 ±17 ADU accumulated in 20 s, 
or 5.95 ± 0.85 ADU s~^ for an object with a Vega-based 
magnitude of j/^'^sa _ 4 35 mag. This fractional flux un- 
certainty corresponds to a systematic uncertainty in the 
Vega-based magnitude zeropoint of 0.15 mag. Using the 
long exposures would not gain much in signal to noise 
ratio, because the SNR is independent of exposure time 
when an object trails in an image. 

The rate of motion of a star at a known right ascension 
and declination allows for the mapping from pixels to arc- 
seconds in the region near the star. The data from Table 
[3] were also used to this end. p Ceti was at a declination 
of —12.29 °, implying an apparent motion across the sky 
of /X = 15" s-i cos 12.29 = 14.66" s'^ We detected the 
object as moving 8.6 pixel in 85 s, or 0.101 pixel s~^. 
This implies a plate scale of 144.9" pixel" ^, so that each 
pixel subtends an area of 20, 995 arcsec^, or 5.8 arcmin'^. 



Source 



Error^ (mag) 



Sky area pixel" ^ 
Vega-system zeropoint 
AB system zeropoint 
Stray & scattered light 
Photometric conditions 
Passband errors 



0.05 
0.15 
0.05 

Unknown 
0.01 

Unknown 



^ Estimated upper limit of systematic 
error. 

We consider this plate scale to be determined with a 
fractional uncertainty of about 5%, due to both the 1 s 
granularity in the timing between exposures and the cen- 
troid measurement uncertainty for the star. Our dy- 
namic plate scale determination is also within 6% of 
the 153.4" pixel" that was returned by astrometry.net, 
based on matching sources across the field to cataloged 
objects. In one of our longer exposures we measured the 
trail length of a star at this same location on the array 
as spanning 31.2 pixels in a 300 s exposure, again con- 
firming the plate scale determination at the 5% level. 

The calibration source, and the region around it where 
we determined the sky brightness values, reside near a 
corner of the field. The rays that are brought to a focus 
at this location pass through the filter at an angle of 
18.6°, while the median LSST incidence angle is 18.9°. 
We therefore consider the center of our effective passband 
at this location to be a good match to LSST's, but the 
shape at the passband edges will still differ somewhat. 

With the photometric calibration and plate scale in 
hand, we were in a position to convert the sky bright- 
ness values, S, into Vega-based units of mag arcsec"^. 
We converted the sky brightness values into units of 
ADU s~^ arcsec"^ by dividing by 20,995 arcsec^ pixel" ^. 
We then computed the Vega-based sky brightness using 
the ffux from the comparison star, as 



-^Vcga 



-2.5 log(0/5.95) -1-4.85 



mag 
arcsec^ 



(2) 



where is the sky flux in ADU s""'^ arcsec"^. Note 
that we selected a region for sky brightness determi- 
nation that surrounded the calibration star to mini- 
mize our sensitivity to large scale fiatfielding errors and 
pixel scale variations. The final zeropoint that converts 
all temporal results of Figure [5] to Vega magnitudes is 
E = AS -I- 17.8 mag arcsec~2 (Vega). 

We also supply the conversion of Vega magnitudes 
through the y band to AB magnitudes, because, while 
the AB and Vega magnitude scales are matched at the 
V band, they monotonically diverge at longer wave- 
lengths. The Spitzcr fiux conversion utilitjQ delivers 
^gVcga _ Q g Ya&g with an uncertainty of 0.05 mag,, 
where y^^ = y^'^sa ^_ AB^^^a^ rp^-^ j^^^j uncertainty 
is adequate for our purposes here. The conversion also 
agrees within the uncertainties wit h the AB-to-V ega off- 
set terms presented in Table 7 of iHewett et al.l (p006). 
Thus the conversion E — AS -I- 18.3 mag (AB) brings all 
temporal results of Figure [5] into the AB system. Table 
[4] presents a summary of the potential sources of system- 

^ At |http: //ssc . spitzer.caltech.edu/wanimiission/propklt/pet/magtojy/ 1 
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atic error that might afflict these results. 



7. CONCLUSIONS 

We have performed first measurements of 7/3 band sky 
brightness in Chile from Cerro Tololo, and have argued 
that these results to apply to the nearby LSST site on 
Cerro Pachon. We observed rms spatial variability of 3- 
4% due to atmospheric gravity waves on scales of a few 
degrees and larger, and a factor of 2 variabihty in flux 
over the course of the nights we observed. These rela- 
tive results make no reference to absolute calibrations, 
and we argue that they apply to nearly any y band vari- 
ant because the sky background is overwhelmingly due 
to narrowband OH emission. Finally, we estimate an 
absolute zeropoint and find the mean sky level is about 
17.8 mag arcsec"^. Our absolute calibration suffers up 
to 20% systematic uncertainty from zeropoint and pixel 
scale measurement errors, while our relative photometry 
is uncertain at about a percent. 

Previous aeronomy work has characterized the OH 
background in depth, and suggests that the variability 
is somewhat predictable. We expect the sky background 
through any similar, CCD-based y band variant to fluc- 
tuate significantly, although nominal variation and ab- 
solute levels will differ depending on the precise nature 
of the total filter response curve. Our qualitative results 
and conclusions should apply to any y band. 

The variability suggests that sky survey projects would 
benefit from a y band sky brightness monitor, and an 
adaptive strategy that would exploit times when the y 
sky is dark to take observations in this band. The back- 
ground in y appears uniform to within 4% over much of 
the sky, so we advise it is better to chase times of low 
background than to find directions of dark sky in y. The 
coherent spatial structures in y band airglow present a 



potential difficulty to sky-based illumination correction 
strategies. Assuming random Gaussian fluctuations, 15 
images or more are needed to achieve sub-percent uni- 
formity over fields of view larger than about a degree. 
The spatial coherence of the gravity waves means the 
fluctuations are not random, making this a lower limit. 

We have demonstrated the effectiveness of this appara- 
tus. Given the observed variability in 7/3 band sky bright- 
ness, there is considerable merit in implementing a long 
term monitor at the LSST site. Having long term sky 
brightness data in hand early will inform the optimal 
scheduling of the LSST system. 
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